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Abstract. The newly introduced Transition Distribution Amplitudes (TDAs) are discussed 
for the 7r-7 transitions. Relations between 7r-7 and 7-7r TDAs for different cases are given. 
Numerical values for the 7r-7 TDAs in different models are compared. GPD's features are 
extended to TDAs and the role of PCAC highlighted. We give hints for the evaluation of 
cross sections for meson pair production in our approach. 
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1. INTRODUCTION 

One of the main open questions in particle physics is the understanding of the 
structure of hadrons in terms of quarks. An important tool for such a purpose is 
provided by hard processes. The large virtuality involved in these processes 
allows the factorization of their amplitudes into hard and soft contributions. The 
hard contribution to the scattering amplitude is known from perturbative QCD 
but the interesting quantity unveiling the structure of hadrons is the soft or non- 
perturbative contribution. In Deep Inelastic Scattering (DIS) this nonperturbative 
quantity is expressed in terms of the parton distribution functions (PDF). PDFs 
can be expressed as the Fourier transform of a bilocal current between the same 
initial and final hadronic state. Generalized parton distributions (GPD) [H, 0] 
extend this concept to off-diagonal matrix elements of the same currents. GPD gov- 
ern exclusive processes with the same hadron in the initial and final state in the 
soft part of the process but with momentum transfer different from zero. Deeply 
virtual Compton scattering (DVCS) is a tipical example of processes governed by 
GPDs. Recently a further "generalization" to transition distribution amplitudes 
(TDA) in which the initial and the final state in the soft part of the amplitude are 
different has been introduced. They have been defined for processes like hadron 
annihilation into two photons and backward VCS in the kinematical regime where 
the virtual photon is highly off-shell but with small momentum transfer t [3] 
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FIGURE 1. Factorization of the amplitude for the process ([T]). 



with H a hadron, or exclusive meson pair production in 7*7 scattering in the same 
kinematical regime [sj 

727^M±7r^ , (2) 

M being either pi or tt. Advocating [3] that the factorization theorems for exclusive 
processes can be extended to the processes under consideration, i.e. ([1]) and ([2]), 
the amplitude can be factorized as shown in Fig. [1] with the TDAs, describing 
the TT — > 7 transition, being the Fourier transform of the matrix element of bilocal 
currents at a light-like distance. 

The nonperturbative nature of the distribution functions imposes the use of 
effective theories, models or phenomenological parametrizations. In Ref. ^ the 
calculation of the 7r-7 TDAs in a covariant Bethe-Salpeter approach has been 
defined. All the invariances of the problem are hence preserved, e.g. gauge and 
translational invariance. For the numerical evaluation, the Nambu - Jona-Lasinio 
(NJL) model for the description of the pion has been used. The Pauli-Villars 
regularization scheme has been chosen because it is Lorentz invariant. 

Estimates for the 7r-7 TDAs have been given in Ref. 01 and a calculation has 
been performed in the Spectral Quark Model (SQM) [8|. Both studies parametrize 
the TDAs by means of double distributions [if. A detailed comparison of the SQM 
and the NJL models in the determination of the pion GPD can be found in Ref. [9|] . 
Finally in Ref. [10] TDAs have been calculated in a non-local chiral quark model, 
confirming the results of the previous calculations. 

In the following Section we introduce the vector and axial TDAs emphasizing 
their connection to the pion transition form factors, appearing in the radiative 
pion decay The presence of an additional contribution due to PCAC is 

explicitly shown. Numerical results and a comparison between the different models 
are presented in Section [3l 

2. DEFINITION OF THE PION-PHOTON TDAS 

General arguments, such as Lorentz invariance, lead to some important properties 
of GPDs. Taking their first Mellin moment, one can relate GPDs to the corre- 



FIGURE 2. Pion pole contribution between the axial current (the cross) and the incoming 
pion-photon vertex. 



spending form factors through the sum rules. Also their higher Mellin moments 
are polynomials in the skewness variable ^ by Lorentz invariance. 

Similarly, as a consequence of Lorentz invariance, TDAs are constrained by 
sum rules and polynomial expansions. The first Mellin moments of n-'j transition 
distribution functions are related to the vector and axial-vector transition form 
factors, Fy and Fa, through the sum rules. The definition of these form factors is 



given from the vector and axial- vector currents [11 
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with = 92.4 MeV, e^^^s ^ i ^ (ri-irs) /2. All the structure of the 

decaying pion is included in the form factors Fy and Fa- The vector current only 
contains a Lorentz structure associated with the Fy form factor. The axial form 
factor Fa also gives the structure of the pion but contains additional terms required 
by electromagnetic gauge invariance. The second term on the right-hand side of 
Eq. (jl]), which corresponds to the axial current for a point-like pion, contains a pion 
pole coming from the pion inner bremsstrahlung: the incoming pion and outgoing 
photon couple with the axial current through a virtual pion (Fig. [2]) as required by 
the Partial Conservation of the Axial Current (PCAC). The third term in Eq. (jlj) 
is a pion-photon- axial current contact term, proportional to fnOfiu- 

Going to the TDAs we introduce the light-cone coordinates = {v^±v^) / 
and the transverse components v-^ = for any 4- vector f^. We define P = 

{p + p') /2 and the momentum transfer, A = p' — p, therefore = m^/2 — 1/4 and 
t = A"^. The skewness variable describes the loss of momentum in the light front 
direction of the incident pion, i.e. ^ = {p — p')^/2P^. Its value ranges between 
t/ {2m\ — t) < ^ < 1. Negative values of the skewness variable can be allowed. 
Regarding the real photon polarization, e'", we have the transverse condition e.p' = 
and an additional gauge fixing condition. We assume that this condition is such 
that e^/P~^ is kinematically higher twist. The standard gauge fixing conditions, 

= or = 0, both satisfy the previous requirement. To leading twist, the TDAs 



are therefore defined 
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witli e(^) equal to 1 for ^ > 0, and equal to —1 for ^ < 0. Here V{x,^,t) and A{x,^,t) 
are respectively the vector and axial TDAs. Hence the axial matrix element contains 
the axial TDA and the pion pole contribution that has been isolated in a model 
independent way 0, @, l2| • The latter term is parametrized by a point-like pion 
propagator multiplied by the distribution amplitude (DA) of an on-shell pion, 
(j){x). Notice that the pion DA obeys the normalization condition J^dx(f){x) = 1; 
the connection through the sum rules of Eq. ([5]) with Eqs. ([3]) and (jl]) is therefore 
obvious. 

The contribution of a pion pole is not a new feature of large-distance distri- 
butions. TDAs like GPDs are low-energy quantities in QCD though their degrees 
of freedom are quarks and gluons. One thus expects chiral symmetry to manifest 
itself, what implies a matching between the degrees of freedom of parton distribu- 
tions and the low-energy degrees of freedom such as pions. Actually, in the region 
X G [— 1^1, 1^1], the emission of a qq pair from the initial state can be assimilated to 
a meson distribution amphtude. 

Here we have defined the TDAs in the particular case of a transition from a vr"*" 
to a photon, parametrizing the processes given by Eq. ([T]). Symmetries relate the 
latter distributions to TDAs involved in other processes. For instance, we could 
wish to study the •y-n^ TDAs entering the factorized amplitude of the process 

Unifying our notations with the notations of Ref . , we define the 7-71^ TDAs 
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FIGURE 3. Diagrams contributing to pion-photon TDA, shown here for an active u-quark. 
The diagram on the right corresponds to a pion "rescattering" . 
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where we have preserved the definition ^ = {p — p')^ / {p + p')^ given before the 
Eq. (p. 

Time reversal relates the n'^-'-f TDAs to 7-71"'" TDAs in the following way 



(7) 



where D = V,A. And CPT relates the presently calculated TDAs to their analog 
for a transition from a photon to a tt^ 



(8) 



3. DISCUSSION 



Basic properties of GPDs and TDAs, like sum rules and polynomiality are related 
to Lorentz and gauge symmetries. Therefore, we need a method of calculation which 
preserves these properties. The main problem here is the description of hadrons, 
as bound states of quarks, preserving these symmetries. One solution is to use 
a field theroretical formalism, with a covariant Bethe-Salpeter approach for the 
description of the hadrons. In this formalism, GPDs and TDAs are integrals over 
the Bethe-Salpeter amplitudes. This method has been developed in Refs. [g], [13] 
for local lagrangians and in Refs. 14, 1^ for non-local lagrangians. In the case 
of pions, the simplest realistic model which realizes these ideas is the NJL model 
within the Pauli-Villars regularization scheme. We will therefore use this model for 
the discussion of the calculation. 

We consider that the process is dominated by the handbag diagram. Each TDA 
has two related contributions , depending on which quark {u or d) of the pion is 
scattered off by the deep virtual photon. The leading contributions to the handbag 
diagram are depicted in Fig. [3] for an active w-quark, with the diagram on the 



right corresponding to a coupling of the bilocal current to a quark-antiquark pair. 
The vector TDA receives contribution only from the first type of diagram, i.e. 
the diagram on the left of Fig. [31 In the case of the axial TDA, a contribution in 
the — < X < \C,\ region arises from the second diagram of Fig. [3l This second 
contribution comes from the re-scattering of a qq pair in the pion channel. It 
contains the pion pole which, according to Eq. ([5]), must be subtracted in order to 
obtain the axial TDA. 

We can express both, V{x,^,t) and A{x,C,,t), TDAs as the sum of the active 
M-quark and the active (i-quark distributions. The first contribution will be pro- 
portional to the d's charge, Qd, and the second contribution to the m's charge, Q„. 
Therefore, we can write 

{x, ^,t) = Qd (x, e, t) + g„ e, , (9) 

with D = V,A and d = v,a. Isospin relates these two contributions. For the vec- 
tor TDA, we obtain = v^^^{—x,^,t). For the axial TDA we have 

aj^-(x,.^,t) = — aj^^^(— where the minus sign is originated in the change in 
helicity produced by the 75 operator. The d'^^^{x,^,t) contribution is non- vanishing 
in the region — 1^| < x < 1 while d'^^-{x,^,t) in the region — 1 < a; < |^|. Given the 

relation (Q, the support of the whole TDA, V"^ {x,^,t) or A^"^ {x,C,,t), is therefore 
X G [—1,1], as required. 

In Fig. m we show the vector and axial TDA calculated in the N JL modeQ for 
different values of the momentum transfer t. The mass of the pion, being either 
= MeV in the chiral limit or = 140 MeV, does not significantly infiuence 
the result. A ^-symmetry is observed for V{x,C,,t) in the chiral limit: the vector 
TDA is an even function of the skewness variable so that we show the results for 
positive ^ only. At the contrary, the shape of the axial TDA depends on the sign of 
the skewness variable. The two distinct behaviours are shown in Fig. |H For positive 
^, the contribution coming from the second diagram of Fig. [3] is dominant and 
produces the maxima around x = 0. For negative values of the skewness variable, 
the contribution of both diagrams have opposite signs, as shown in Fig. [5l Given 
Eq. iQ, we can say that isospin relates the value of the vector and axial TDAs in 
the \^\ < X <1 and — 1 < a; < —\^\ regions, 

V{x,C,t) = -^V{-x,^,t) & A{x,^,t) = ^A{-x,^,t), \^\<x<l . (10) 

The factor 1/2 corresponds to the ratio between the charge of the u and d quarks. 
We observe in Fig.lHthat our TDAs satisfy these relations. It must be realized that 
the relation (ITO!) cannot be changed by evolution. 



For numerical predictions, the standard values of the parameters given in Ref. [TB| are used. 



Vector TDA (m^=140 MeV and t=-0.5 GeV ) 



Axial TDA (m^=0 MeV and t=-0.5 GeV ) 




FIGURE 4. Vector and axial TDA in the NJL model. The amplitudes are lower for higher 
{—t) values, as it can be inferred from the decreasing of the form factors with {—t), connected to 
the TDAs through the sum rules. 



The obtained TDAs should obey the sum rules, as already mentioned. For the 
form factors given in the NJL model (see Appendix of Ref. [6[]), the sum rules 

f dxD^\x,^,t) = ^Fn{t) , (11) 

with D = V,A, are recovered. In particular we obtain the value Fy^ (0) = 0.0242 
for the vector form factor at t = 0, which is in agreement with the experimental 
value Fv{0) = 0.017± 0.008 given in [17]. We also obtain F^^(O) = 0.0239 for the 
axial form factor at t = 0, which is about twice the value FJ^ (0) = 0.0115 ±0.0005 
given by the PDG [12]. 

We expect the TDAs to respect the polynomiality condition. However no con- 
straint from time reversal enforces the polynomials to be even functions of the 
skewness variable, i.e. the polynomials include all powers of ^ 

where C^^{t) are the generalized form factors. The relation (1121) has been numer- 
ically verified. In the chiral limit, we have numerically found that the odd powers 
in ^ go to zero for the polynomial expansion of the vector TDA. 

Other studies of pion-photon TDAs have already been done 0, 0, 
[sl , double distributions have been used. Therefore polynomiality is implemented by 
definition and cannot be considered a result. The aim of the author of Ref. 0] is to 
provide some estimates of the vector and axial TDAs on the basis of the positivity 
bounds. In this way we must compare only the order of magnitude of the obtained 
amplitudes, which is similar to ours. 

The vector and axial TDAs calculated in the SQM [8], NJL model [6| and non- 
local chiral quark model [lo!] are compared in Fig. [HI The authors of the first 
reference use the so-called asymmetric notation. The comparison is here awkward 
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FIGURE 5. Contributions to the axial TDA for both positive — 0.25, soHd line) and negative 
= —0.5, dashed line) values of the skewness variable. In each case, and in the x G [— ICU^I] 
region, the contribution coming from the first diagram of Fig. [3] is represented by the dashed- 
dotted lines and the non-resonant part of the second diagram of Fig. [3] is represented by the 
dotted hnes. 



since the authors define, after their Eq. (3), C = {p-i—V-K^^ Ivt while, after their 
Eq. (4), the definition —C, = {p^—p-„).n is given. We nevertheless decide to use 
the standard relation between their asymmetric notations and our symmetric 
ones Their functions V^qm and A^qm corresponds to our d"*^'^ given in Eq. ([9]). 
The normalization condition is different from the one used here and the results 
quoted by these authors must be corrected, for the vector TDA, by a factor 
487r^ ^/2fT^Fv{0)/m-,r ~ 10 before comparison. From Fig. [6]we conclude that there is 
a qualitative and quantitative agreement, for the vector TDA, between the results 
of Ref. and those obtained in the NJL model. Regarding the axial TDA we 
observe, in addition to the normalization factor 487r^ V^/tt -^yi(0)/^«^ ~ 10, a change 
in the global sign due to different definitions. In the first version of the axial 
TDA for positive values of ^ is giverH. It coincides with the results in the NJL 
model, as observed in Fig. [61 Surprisingly, the result presented in Ref. [8| coincides 
with our result for negative ^ (perhaps due to the change in the definition of ^ 
mentioned above). 

In Ref. the TDAs are calculated in three different models. The first one is 
a local model which pole structure has some similarity with the one of the NJL 
model. The two other models, i.e. semi-local and full non-local, follow the results 
of the local one. The most prominent difference between the results obtained in 



There is a typographic error in Eq. (23) of the first version of this reference, where a factor 
My/ 6 must be dropped. 



the non-local and the NJL models is the appearence of important odd powers in ^ 
in the polynomial expansion of the vector TDA. We know from Ref. that, for 
non-local models, there are additional contributions to those calculated in Ref. [lO| • 
In the case of PDFs, disregarding these contributions can produce small isospin 



violations [l5j. It can therefore be considered that, on that point, the results of 



Ref. must be confirmed. For numerical comparison, the results obtained in 



10| must be corrected by a factor 27r^\/2/^Fy(0)/m7r ~ 0.45 due to the use of 



a different normalization condition. For the axial TDA, there is a change in the 
definition of the skewness variable between the caption of Figs. 3, 5 and Fig. 9 of 
Ref. [10]. We observe, see Fig. [6l that our results in the NJL model coincide with 
the results of the latter reference if the convention of the caption of their Fig. 9 is 
chosen. 

It can eventually be concluded that there is no disagreement between the 
different studies concerning the 7r-7 TDAs besides the ambiguity in the definition of 
the skewness variable. Calculations of Refs. j^, l3] are performed in the chiral limit, 
where ^ runs from —1 to 1. The symmetric nature of this interval makes difficult to 
check the sign of ^. On the other hand, the NJL calculation is given for physical 
pion mass. In this case, the kinematics of the process imposes t/ {2m'^ — t) < ^ < 1. 
From Eqs. (22) and (23) of Ref. we observe that there is a pole in the axial TDA 
for the limit value C, = t/{2m^ ~t), preventing us from going through unphysical 
values of ^. Moreover, the sum rule Eq. ( ITTi) for both the vector and axial TDAs 
are here satisfied for physical values of ^, and broken in the unphysical regions 
^ < t/{2m^ — t) and > 1. We therefore conclude that the choice of sign in Ref. [sj 
is consistent and gives a guideline for comparing with other models. 

The 7r-7 vector and axial Transition Distribution Amplitudes have been defined. 
The results in different models have been discussed. In particular, the numerical 
results of Ref. [8|, i6j, ilOj] have been compared. The use of a fully covariant and gauge 
invariant approach guaranties that all the fundamental properties of the TDAs will 
be recovered, in particular, the right support, i.e. x G [—1,1], the sum rules and 
the polynomiality expansions. Hence, in the NJL model, these three properties are 
not inputs, but rather results of the calculation. 

Recently, cross section estimates for the process ([2]) have been proposed in Ref. [sj 
using, for the non-perturbative part, t-independent double distributions, in a first 
approach, and, in a second, the t-dependent results of Ref. 0]. In the fine of 
sight of this reference, similar estimates for the meson pair production in hard 
7*7 scattering using the results for the TDAs cited above could be given [l^. In 
particular, it would be worth investigating the model independence of the structure 
independent terms [l^ and then estimating the pion pole contribution to the cross 
section pJi]. 
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FIGURE 6. Comparison of the 7r+ ^ 7 TDAs for t = -0.1 GeV^ of Refs. % |Tq for = 
MeV and Ref. [fi] for ra^ = 140 MeV. On the left, we have: the vector TDA for ^ = ±0.5 as 
a single (solid) curve for the results of both the NJL model and SQM (these four curves are 
indistinguishable); the result of the non-local xQM calculation for ^ = 0.5 (dotted line) and for 
^ = —0.5 (dashed-dotted line). On the right, we have the results of the axial TDA with the choice 
for the sign of ^ as discussed in the text. 
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